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1 Introduction 

The path integral approach has been extensively applied to calculate the thermodynamic prop- 
erties of the quantum field theories |jl|, |^, |^. Through this approach, it has been calculated the 
leading contribution, in the high temperature limit (/3 <^ 1, where /3 = ■^), to the effective 
potencial of these theories^, Q. We also have a standard high temperature perturbation theory 
derived from the path integral expression of the partition function]^, ^. In the calculation of 
the Helmholtz free energy, via the path integral approach, we do not keep track of its contri- 
bution that are /3-independent. In general this contribution is irrelevant to the thermodynamic 
properties of the model, but even we are not going to deal here with supersymmetric models, 
we remember that this contributions are crucial to verify if the supersymmtry is broken or not 
at finite temperature0, |, ||, since for a theory to be supersymmetric the value of its vacuum 
energy has to be zero. 

For theories involving bosons, in the path integral we first integrate over the conjugate 
momenta. Bernard [^ showed that in this case we have to be carefuU in constructing the path 
integral to get the overall /3-dependent constant of the grand canonical partition function. 
This constant does not appear in pure fermionic models since in these cases the momentum 
conjugate to the fermionic fields are their own hermitian conjugate. Even for these fermionic 
models, along the calculation of the partition function done via the path integral approach, we 
get /3-independent terms that are dropped out (see the free fermion case that is fully discussed 
in reference [^). 

The partition function of free fermions models is calculated directly from its path integral 
expression since it is equal to the determinant of its dynamical operator [^, |l^. When the 
lagrangean density of model has self-interacting fermionic terms, it is not possible any more 
to calculate exactly the non-gaussian non-commutative path integral. One common approach 
is to get the bosonized version of the model|jll|. Another standard way to handle the path 
integral over non-commutative function is to do a perturbation theory taking care of the signs 
coming from fermion loops |jl|, ^, ^. 

The interesting properties of the non-commutative Grassmann algebra has been applied 
to get the contributions from spin configurations to the partition function of the classical 
bidimensional Ising model[0]. In reference |T^ Charret et al. proposed a new way to calculate 
the coefficients of the high temperature expansion of the grand canonical partition function 
of self-interacting fermionic models in (i-dimensions {d > 1). They applied the method to 



the Hatsugay-Kohmoto model[0], that is an exactly solvable model. The approach was also 



applied to the unidimensional Generalized Hubbard model to get the coefficents up to order 



j3^ of the high temperature expansion of its grand canonical partition function ||14||. Differently 
from other approaches, these coefficients calculated by this method are analytical and exact. 

Up to now the approach of Charret et al. has only been applied to fermionic models already 
regularized on a lattice with space unit one|TB|, |T1|. The aim of the present paper is to use this 



method to calculate the Helmholtz free energy of the free fermion Dirac, that is a continous 
theory and whose exact result is already known in the literature 0. It is also important to 



check its application to a fermionic lattice model ||T0| 



In section 2 we give a summary of the results of the method of Charret et al. fully described 



in reference ||13[. In section 3 we apply the approach to the free Dirac fermion using two 
expansions: in subsection 3.1 we first expand the fermionic field operators in the basis of 
the eigenstates of energy and in subsection 3.2 we consider the naive fermion model on a 
three dimensional space lattice. In both ways, we show that is possible to re-sum the high 
temperature expansion of the grand canonical partition function of the model and compare our 



results with the known ones in the literature^, |T0|- In appendix A we present the formulae 
related to the lagrangean density of the free Dirac fermion. In appendix B we have a dictionary 
of some formulae showing their continuun and discrete expressions. 

2 A Survey of the Approach of Charret et al. 

The grand canonical partition function of any quantum system can be written as a trace over 



all quantum physical states P, |T5 



Z(/5)=Tr[e-^K], (1) 

where P = -^, k is the Boltzmann's constant and T is the absolute temperature. The operator 
K is defined as K = H — /iN, H being the hamiltonian of the system, fi is the chemical potential 
and N is a conserved operator. In the high temperature limit (/5 <^ 1), Z{j3) has the expansion 

oo (_iQ)n 

Z{(3) = Tr[l] + Y.^-^Tr[K-]. (2) 

The spectation value of any fermionic operator can be written as a multivariable integral over 
Grassmann variables p!6|, 0. The mapping among the fermionic operators that appear in the 



fermionic model and the grassmannian generators is such that they satisfy the same algebra. 
Let aj and a^ be the hermitian conjugate fermionic operators that satisfy the anti-commuting 
relations 



{ai,a]} = 16ij and {a.i,a.j} = 0, (3) 

where i = 1,..,^/. The generators of the associated Grassmann algebra has dimension 2^-^'^, 
and can be written exphcitly as {%,■'■ ;Wj '71; '■' ;W}- They satisfy the following anti- 
commutation relations: 



The mapping 



{rii, rij} = 0, {rii, rjj} = and {rn, rij} = 0. 



Vi 



and 



_d_ 



(4) 



(5) 



preserves the algebra (H) due to the algebra (§) satisfied by the grassmannian generators. 

For any self- interacting fermionic model in a (i-dimension lattice {d > 1), the coefficients of 
high temperature expansion (Q) can be written as the multivariable Grassmann integral |]I3 



Tr[K"] 






Y[ dT]idr]j e^ 

x/C®(r/, 7],u = 0)/C®(r7, r],u=l)... }C®{r], r],u = n-l), 



1=1 



(6) 



A^"^ is the number of sites in the d- dimensional lattice. Matrix A is independent of the operators 
H and N and is equal to 



A 



ATT Q 
Q A^T 



(7) 



Each element of matrix A is a matrix of dimension nN'^ x nN'^ and Q) is the null matrix in 
this dimension. The indices /, J in matrix A vary in the interval, I, J = 1,2,..., 2nN'^. The 
matrices A^T and A^T are equal and 



A 



TT 



U 



Q)j^dy.j^d l^YdxATd —Ij^d^j^d 



Q)j^d^j^d 



V '^N'ixN'i ^'ixN'i ^'ixN'i ■ ■ ■ '^N'ixNd J 

where Ij^d^j^d and (J)j^d-^j^d are the identity and null matrices of dimension A^*^ x A^'' respectively. 
Due to the fact that the submatrices A^T and A^-T are null, the result of the multivariable integral 
(I) is equal to the product of the independent contributions of the sectors aa =tT and aa =ii. 



The grassmannian function /C®(f/, rj) in eq. (^) is obtained from the normal ordered operator 
K by doing the naive mapping: aj^. — > f// and aj.o- -^ r]i, where i = 1, 2, ..., A^*^ and o" =1, |. 

Finally, the result of the integrals that contribute to Tr[K"] are independent of the aa 
sector since A^^^^ = A^^^. We present here only the results of the integrals of the sector a a =tT- 
All those integrals are of the following type: 



.nN" „_„wd 



n dr],dr]i VhVk, • --VuVk^ e^-^- "'^^ ''^ (9) 

where L = {li, . . . , l^} and K = {ki, . . . , km}- We should remember that the grassmannians 
functions /C® are polynomials of the generators of the algebra. The products rjf] are ordered 
such that li < I2 < ... < Im and ki < k2 < ... < km. The results of the integrals of type (|^) 



are 18 



M{L, K) = (-l)S™ i('-+'=«)a(L, K), (10) 

where A{L, K) is the determinant obtained from matrix A^^ by the cut of the lines /i, I2, ...,lm 
and the columns ki,k2, ...,km. M{L,K) is a co-factor of matrix A^^^^. Result (|TD|) reduce the 
calculation of the multivariable integrals over anti-commuting variables to the calculation of 
co-factors of a well defined matrix whose elements are commuting numbers. 

In general it should not be easy to calculate the co-factors of matrix A°"'^ since each of its 
elements is a matrix of dimension A^'^ x A^*^. However, due to the block structure of matrix A°"°" 



Charret et. al showed in reference |13] that is possible diagonalize it for any value of n and N . 



The results are analytical and allow us to take the thermodynamic limit. 

3 Grand Canonical Partition Function for Free 
Relativist ic Fermions 

The lagrangean density of the free Dirac fermion is: 

£ = ^(f,t)(,9 -m)^(f,t), (11) 

whose fermionic field operators satisfy the anti-commutation ralations 

{*^(f,t),n,(f,t)} = ^5^,5(f-f'), /i,z/ = l,2,3,4 (12) 

and 

{^^(x, t), ^,(f , t)} = {n^(f, t), n,(x', t)} = 0, (13) 



where Iliy{x,t) is the canonical momentum of \E'j,(x, t) and Ili^{x,t) = 2\1/J,(x, t). 

From the lagrangean density (|TT|) and the equation of motion satisfied by the fermionic field 
operators, the total hamiltonian operator of the system can be written as 

H= / rf3xi^Vf,t)9o*(x,t). (14) 

Jv 

We use the natural units where h = c = e = 1. Our metric is diag{g^y) = (1, —1, —1, —1). 

3.1 Expansion in the Basis of Eigenstates of Energy 

The fermionic operator \l/(x, t) written in the basis of eigenstates of energy of the free Dirac 
fermion is 

^(^'^) = 77^EE[ar.(^)«.(^)e-^'=--^ + bt(^)t;,(^)e^^'^-a (15) 

k '■=1 



where k^x^ = k^xo — k.x, and ko = ylk]"^ + iv? > 0. The destruction fermionic operators a^, 
hj. and their respective hermitian conjugate satisfy the anti-commutation relations 

{sLr{k),al{k')} = 5rs5j:p and {b^(fc), bJ(A^')} = (5^,(5^ g,. (16) 

All others anti-commutation relations of these operators are null. The spinors components 



Ur{k) and Vr{k) are given by eqs. (|49D and (M)- 

The hamiltonian operator written in terms of creation and destruction operators of defined 
energy becomes 

H = -1EE^0 + EE koH{k)Bir{k) + hl{k)hr{k)], (17) 

k ^=1 k '■=1 

1 been the identity operator. We define Eq as the vacuum energy, Eq = (0|H|0) = —^J^k^O: 
|0) being the vacuum state of the fermionic model. 

Our aim is to calculate the grand canonical partition function of the free Dirac fermion in 
contact with a reservoir of heat and electric charge. Eqs. (^ , (|^) and ( PUJ ) allow us to calculate 



the coefficients of the high temperature expansion of this function. For the present case, the 
operator K in expression (|l]) is: K =: H : +Eq1 — yuQ, where /i is the chemical potential and 
Q is the total electric charge operator of the free relativistic fermions. Operator Q written in 
terms of the creation and destruction operators is 



Q = 2 5:il + ^^[at(^)a,(^)-bt(^)b,(A^)]. (18) 

k k '■=1 

= gol+:Q: ■ 

Finally, we have that K operator can be written as 

2 2 

K = EE(^o-/u)at(^)a,(^) + ^^(A;o + /i)bt(^)b,(^) + ]l(Eo-/iQo) 
k '■=1 k '^=1 

= Ka + K, + l{Eo-fiQo). (19) 

Due to the anti-commutation relations ([T6|) , we have that [K^, K?,] = 0, and therefore Z{(3) can 
be written as 



where 



Z{f3) = e-^(^°-'^«°)Z,(/3).Z,(/5), 



ZM = Traie-f"^^] and Z,{(3) = Tn[e- 



-/3Kb 



(20) 



(21) 



The calculation of the functions Za{P) and Zi,{[3) are equivalent. We present here only the 
details of the calculation of Za{(3). For a free Dirac fermion, the traces for r = 1 and r = 2 are 
equal. Then, 



= IT Tr^a^ [g-/3(fco-M)ni'> (fc) j 



(1), 
2 



f22) 



In Tr^^ we have the vector k fixed. Since n^^^ (A;) is a commuting operator, we apply the Newton 
multinomial to write 



,(1) 



-py 



Tr(^)[,-mo-.W. '(.)] = Tr«[|] + J2 '-^ih - /.)"Tr«[n«(fcl"]. 



n=l 



n\ 



(23) 



We should note that for fixed k the anti-commutation relations (|T6|) are identitical to the 
relations (^, therefore, Tr'^^[n^^\k)'^] can be written as a Grassmann multivariable integral 
(^) with an equivalent mapping to eq. (^ for the associated generators of the non- commutative 
algebra. The traces that contribute to eq.(p3|) are written as the following anti-commuting 
integrals: 



7=1 

and 



Tr^\^] =jf[ %(/)%(/)eS^-'^^^(')^"''''^(^\ (24) 



-' 1=1 

x%(l)^fc(l)---%H^fcW' for (n>0). (25) 

The grassmannian expression for Tr^^)[n^^^(fc)"] is obtained from eq.(^ with A^ = 1 and d = 1. 
Therefore the elements of matrix A*^^^^ in eq.(|]) are just numbers. From eqs.(^ and (p!0|), 



we conclude that having the 77's in the integrand corresponds to deleting the first n lines in 
matrix A*^^^^ and, in the same way, having the r/'s corresponds to deleting the first n columns 
of the same matrix. Therefore, from the expression of matrix A*^^^-* (see eq.(^), for arbitrary 
n, we realize that the matrix A^^^^ obtained after the cuts of n-first lines and n-first collumns 
is an upper triangular matrix whose determinant is equal to 1. Beside this, we have that 
det[A(^^)] = 2, for any value of n. In summary, for arbitrary n we have that 

rrW[l] = 2 and Tr«[n«(^)1 = 1- (26) 

Substituting results ( pB]) in expression (^) and resumming it we get 



y^(i)[g-/3Ki^)] = 11(1 + e-^('=°-'^)), (27) 

k 

that returning to eq.(PD|) gives us 

Z{(3) = [[]e2^('="+'^ni + e-^('=o-'^))(l + e-^(^«+^))]2. (28) 

k 

The relation between the Helmholtz free energy and the grand canonical partition function 



IS 



ra 



W{(3) = -^Hzm. (29) 



Substituting result (128|) in eq. (|29|) the expression derived for the Helmholtz free energy is 



W{p) = -| ^[/3(A;o + /x) + ln(l + e'^^'^o-^)) + ln(l + e-'5('=»+'^))]. (30) 

^ k 



3.2 Free Dirac Fermions on the Lattice 

The lattice calculation of models including fermions has been an important tool to learn the 
properties of these models. The aim of this subsection is to show that the method of Charret 
et al. can be equally well applied to the lattice version of fermionic models. To do so, we 



consider the most naive lattice realization of the free Dirac fermion|[TO[]. We remember that the 
crucial point to apply the method of Charret et al. is to work with fermionic field operators 
that satisfy the anti-commutation relations (^). 

The hamiltonian operator of the free Dirac fermion is 

11= I Sx z^^(f,t)[7°7-V-m7°]^(x,t), (31) 

J Y 

where 7 = (7^,7^,7'^). From appendix B, we get that the operator H written on the lattice 
becomes 

N 

- 4 3^2 

H = E E EVy*a(^«''^)(7V)a/3[^/3(na + ia,t)-^;3(na-ia,t)]- 

- ia^m^aiina, t){-f^)ap^fs{na, t)} , (32) 

where a is the distance between the nearest sites in each space direction. The space point x on 
the lattice is written as: x = na. N is the total number of space sites in each direction. For 
simplicity we take to be N even. 

Operator K in expression (|1]) is: K = H — yuQ, where ^ is the chemical potential and Q is 
the total electric charge operator. The discrete expression of Q is 

N_ 
2 

Q = a3 ^ ^],(na,t)^^(na,t). (33) 

ni,n2,n3=-^- 

Imposing the space periodic boundary condition for the fermionic field operators, their 
Fourier decompositions are 

vl/„(na,t) = ^ Y: ^„(!l^,t)e^f^^-, (34) 

where a = 1,2,3,4 and L = ^. The fermionic field operators \1'q,(x, t) and \E'J^(x, t) satisfy 
the anti-commutation relations ([T^) and (13) which imply that the Fourier components of the 



fermionic field operators obey the relations: 



,7r 



,iTk 



.Til 



,7ck 



{MY,t)Jl{^,t)} = 5afs 6'"! and {Mf,t),M^,t)} = 0, (35) 

that are identical to relations @, the necessary algebra to apply the method of Charret et ai. 
In momentum space, operator K is written as 



K 



1 
a 



h,l2,h = -^ 



where 



and 



K,.^t(p:,)Kgr).(^M, 



(36) 



(37) 



The matrix R ( ^/j is defined as 



^(|z:t 






-'1^ 




(38) 



(39) 



- yu)ll sin (jlja) <Jj 
ijlja) cTj —a{m + yu)! 

We have an implicit sum over j in the off-diagonal elements of matrix R and 1 is the identity 
matrix of dimension 2x2. Due to the anti-commutation relations ([35| ) we have that [K^ K^] = 0. 
Therefore, the grand canonical partition function of the model becomes 



zm 



— -1 

2 '- 

n 



Trr 



e ^ ' 



(40) 



where Tr^ means that the trace is calculated for fixed /. We should notice that the operator 
Kj-is not diagonal in momentum space. We make the similarity transformation: PRP~^ = D, 
where the diagonal matrix D is 



D 



(D A_ll 



(41) 



where A+ and A_ are the eigenvalues of matrix R, and 



-a{fi ± ui{l)) 



(42) 



10 



with 



u{l 



;r 2 I ~ 2 



^ 2 



m^ + pi + p2 + P3 



(43) 



and / = {h-ih-ih)- We followed reference [10 and made the change of variables 



Pi = -sin{pia) 
a 



(44) 



where pi = jU. The new fermionic field operators \I'' = P\I' and \E' = \E''''P ^ also preserves 



the anti-commutation relations (PH). 

The function Z{f3) written in terms of the new fermionic field operators has the same form 



as the r.h.s. of eq. (|22|) . Following similar steps, we get the Helmholtz free energy for the model 



on the lattice, that is. 



W(/3) 



-2 Y. (f^+^in 



h,h,h=— — 



_1 ^ [/^(^i + e-^(^(^)+^))+/n(l + 

Pill— N 



-/3(<i(0-/^) 



(45) 



In the limit a ^ the function VV(/3) agrees with the result of reference ||T0| and is equal to 
twice the result (|30D. Here as well in the usual lattice calculation with fermions, the doubling 
problem is lifted by including the Wilson term 



2 Jv 



(46) 



in hamiltonian (|3T| ) and r is the Wilson's constant. 



4 Conclusions 



Recently Charret et al. proposed a new way to calculate the coefficients of the high temperature 
expansion of the grand canonical partition function Z{(3) of any self-interacting fermionic model 
in c?-dimension {d > l)|l^. To apply this method to calculate these coefficients, it is enough to 
write the second quantization expression of the hamiltonian and a conserved operator in terms 
of operators that satisfy the anti-commutation relations (^. In this approach, at each order (3^ 
of the high temperature expansion of the function Z{f3), the calculation of the coefficients is 
reduced to get the co-factors of a matrix with commuting entries (see matrix A in eqs. (|^) and 



11 



(^)). For a fixed number of sites on tlie lattice, all the mathematical objects in the calculation 
are well-defined. This approach avoids to calculate the fermionic path integral of Z{(3). 

We applied the method of Charret et al. to the free Dirac fermion by first expanding it in 
the basis of the eigenstates of energy of the free fermions. The hamiltonian and the total electric 
charge operators (eqs.(0) and ([T8|) respectively) include the respective vacuum contribution. 
Result (|30|) gives two divergent terms to the Helmholtz free energy W(/?): the vacuum energy 
and the electric charge of the vacuum. This last divergent term does not appear when we do 
the calculation via the usual path integral approach^, while the /3-dependent terms in W{(3) 
are identical in both methods. In the path integral calculation divergent terms that are j3 
and /i-independent are dropped out. This does not happens in the present approach. 

Since the calculation of fermionic lattice models are an important tool, we also applied the 
approach of Charret et al. to a naive version of the free Dirac fermion on the lattice [|10|- In 



this case we also get the doubling problem that is lifted by the Wilson term, as happens in the 
usual fermion lattice calculations. 

In summary, we can affirm that the method of Charret ei al. can also be applied to continuous 
fermionic models. It is an analytical approach that allows without ambiguity to calculate all 
terms of the function Z{I3) for the free Dirac fermion, including the divergent terms coming 
from the vacuum contribution. The next step is to study in this approach the renormalization 
scheme associated to physical quantities for fermionic models with self-interacting terms. 
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Appendix 

A Useful Formulae 

In the lagrangean density for the free Dirac fermion (see eq.fpJ])), 

£ = 'l'(f,t)(^ -m)^(f,t), (47) 

the Dirac matrices 7^, /i = 0, .., 3, are: 

12 



(4J 



The matrices cxj, i = 1,2,3 are the Pauh matrices, and 1 and 0are the identity and null matrices 
of dimension 2x2 respectively. We have ^{x,t) = — ■i\l/J,(a;, t)7o. 
The spinors components Ur{k) and Vr{k) in eq.(|T3]) are: 

/ 1 \ / \ 





Ui{k) 



Im + ko 
2kn 



fc3 



m+ko 

. ki+ik2 . 

\ m+ko / 



U2{k) 



Im + ko 
2/Cn 



1 

ki—ik2 

m+ko 

I -fc3 , 

\ m+ko / 



(49) 



t;i(fc) 



/m + /Cn 



2A;n 



/ fel-ife2 


\ 


m+ko 
m,+ko 



V 1 


/ 



V2{k) 



Im + ko 



2kn 



( -^ ^ 

m+ko 

ki+ik2 

m+ko 

1 





(50) 



B Summary of Expressions on the Lattice 

We assume that each direction in the three dimensional space has A^ sites. For simplicity A^ 
is assumed to be even. The distance between nearest points along each direction is a. In the 
space lattice, each point is written as a; = na, where n = {ni,n2, n^), rii = —^, —^ + 1 



2 ' 



' ' ' ' ' 2 ' 



and i = 1,2, 3. 

We have the following dictionary between the continous to the discrete representations of 
the formulae: 



d^ X 



V 



6{x — x') 



JV 
2 

n\,n2,nz=-^- 

■},^n,rh — o'-'?ii,m,i'-'n2,m2"n3,'"3' 



where x = na, x' = ma with n = {ni,n2, n^) and m = (mi, m2, ms) respectively. 



Finally, the first derivative defined on the lattice is defined as[]TO 



(51) 
(52) 



(9j\I^Q,(s, t) ^ — [^a{na + ia, t) — ^^.(na — ia, t)] , 



where i is the unit vector in the direction i. 



(53) 
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